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METHOD OF REPRESENTATION OF ACOUSTIC SPECTRA AND REFLECTION 
CORRECTIONS APPLIED TO EXTERNALLY BLOWN FLAP NOISE 

by Jeffrey H. Miles 
Lewis Research Center 

SUMMARY 

A computer method for obtaining a rational function representation of an acoustic 
spectrum and for correcting reflection effects is introduced. The reflection effects 
were estimated from one -third -octave sound pressure level (SPL) data. The reflection 
and atmospheric attenuation effects were taken into account at the time the power spec- 
tral density function needed to represent the measured SPL wa^ selected. The proce- 
dure is based on minimizing the square of the error between the measured SPL and the 
calculated SPL. 

The most important results obtained were that (1) a rational polynomial with a first- 
order numerator and a third -order denominator was a satisfactorily expressed Strouhal 
normalized mean -square pressure spectral density level and (2) the N independent 
source model of P. Thomas was satisfactory in calculating the reflection effects when 
this model was extended by the use of a reflected -ray transfer function. The computer 
procedure was applied to noise from a full-scale externally blown -flap system having a 
three -flap wing section with a quiet 6:1 bypass ratio turbofan engine having a coannular 
nozzle. The three -flap wing section simulates the take-off condition of a conceptual 
STOL aircraft. 


INTRODUCTION 

Short takeoff and landing (STOL) aircraft may become an important means of trans- 
portation in urban areas. But, strong adverse public reaction to aircraft noise means 
that future STOL aircraft may be required to meet Stringent noise level regulations. 
Thus, noise will be one of the primary factors in the selection of a STOL propulsion- 
lift system. One candidate source of propulsive- lift is the externally blown- flap (EBF) 
system. 



Each EBF STOL concept introduces noise sources not present in conventional take- 
off and landing (CTOL) airplanes. These EBF STOL propulsive- lift concepts are being 
investigated by both scale -model and full-scale engine tests. Concurrently, theoretical 
and experimental studies are being conducted on various noise source mechanisms, and 
empirical prediction methods are being devised. 

In a recent study conducted at Lewis (ref. 1), a rational function representation of 
flap noise spectra was shown to be useful (1) as an empirical equation to fit EBF STOL 
spectral data, (2) as a function that can be related to certain noise- source mechanisms, 
and (3) as part of a computer method for making corrections to acoustic reflections in 
static test data. An objective of this report is to provide the details in the development 
of the equations^ described in reference 1. 

In spite of many advances in the understanding of noise sources, procedures for 
calculating noise from EBF STOL airplanes starting with basic aeroacoustic mechanisms 
are not currently available. Hence present design studies depend on model and full- 
scale tests and on empirical relations determined from these tests. 

For design studies a most important aspect of the measured noise is the one -third - 
octave sound- pressure level spectrum designated in this report by SPL. It is used to 
compute the perceived noise level (PNL) and the effective perceived noise level (EPNL), 
which are the criteria presently used to judge the effect of aircraft noise on the com- 
munity. 

The empirical rational function discussed in reference 1 may be used to represent 
broad-band SPL spectrum data. With the representation described, sound spectra can 
be stored in a computer in a compact form and thus facilitate comparison of large 
amounts of data from different experiments. The representation can also be used to 
extrapolate data to untested conditions and is useful in making estimates of PNL, EPNL, 
and footprint calculations . 

Another objective of this report is to provide an example of the application of this 
empirical rational function representation to a set of data. As an example, the proce- 
dure is applied to data from a full-scale externally blown- flap (EBF) system with a 
three -flap wing section. The EBF propulsion system is operated to simulate the takeoff 
condition. The wing section is positioned to simulate the takeoff condition of a concep- 
tual STOL airplane (ref. 2). The data are from an acoustic program conducted under 
the direction of Lewis using a quieted 6:1 bypass ratio turbofan engine having a coannular 
nozzle (refs. 3 and 4). Reflection effects evident in this research program make it a 
good example to demonstrate the use of the procedure. 
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ANALYSIS 


APPROACH 

The approach taken to represent the measured flap noise spectra is based on a 
model that assumes the measured one-third-octave band SPL spectrum can be repre- 
sented as the sum of a lossless spectrum SPL'*’, a ground reflectance | t | and an 
atmospheric attenuation level Aj. Thus, SPL = SPL'*' + 10 logjQ|T|^ + A.^,, dB refer- 
enced to 2x10 ^ Pa^ (2X10”^ u bar). The same approach was used in reference 1 where 
a condensed description appears. 

The approach taken to represent the lossless spectra will be discussed first. Since 
the development is long, references will be made to appendixes for the details. Next, 
the equations necessary to relate the lossless SPL to the SPL with ground reflection 
effects will be discussed. Finally these results will be used in presenting the complete 
model. 


LOSSLESS SPECTRAL DENSITY 

A semiempirical method of representing SPL spectra was briefly described in ref- 
erence 1 and is presented herein in greater detail. Other procedures for describing 
externally blown- flap SPL spectra are described in references 5 and 6. But, these 
methods are not as universally applicable as the method described. 


Spectral Density Function 


The approach taken to represent the lossless flap noise spectra is based on finding 
a function of Strouhal number such that the one -third-octave lossless sound pres- 

sure level (SPL"**) can be fitted by 


( 1 ) 


where 


f 

r 1 


(IoOASPL/10) 

D^ 

Afc(n— 


1 

L 



OASPL = 10 log Yj 
I 

1 2 
■^1 Pascal (symbol Pa) equals 1 N/m . 


( 2 ) 


3 



and where is an effective diameter (ref. 7), is an effective velocity (ref. 7), 

is the band width of the one -third -octave band having a center frequency f^,, and 
I is the one -third -octave frequency band. (All symbols are defined in appendix A, ) 
Specifically, 


D. 

^ = j27T (f^ — 1= j27T St 
e> 


(3) 


where 


and St is the Strouhal number. 

The semiempirical method to be discussed uses a certain function y?(n) to repre- 
sent the lossless flap noise spectra. The function was developed to represent spectra 
having only a single peak. Thus, for example, the function would be applicable to repre- 
sent single and coannular nozzle jet noise data and EBF noise produced by a single or a 
coannular jet blowing on a flap. Other requirements on the function follow: 

(1) The function should be tailored for applicability to large Strouhal numbers. 

Curve fits that are accurate at large Strouhal numbers are necessary, since EBF en- 
gines are likely to have large effective diameters and low velocities. 

(2) It should be consistent with theoretical considerations. 

(3) It should depend on as few mathematical coefficients and operations as possible 
since this will reduce computer calculation time. 

(4) It should be applicable to a wide range of system configurations (e.g. , single or 
coannular nozzles, over-the-wing: and nnder-the-wing). 

From equation (1) 10 log|^Q0^(J2j| becomes 


lOlog^Q [^(S^)] 



D„ 

f2=j27rf^(0 ^ - 10 log^Q 

Af„(n — 

L 


(4a) 


The quantity represented by equation (4a) is the Strouhal normalized mean-square pres- 
sure spectral density level, hi this report, expression (4a) will be referred to as the 
PSD'*’ function of Strouhal number. Specifically, PSD'*’ is calculated from or represents 
SPL data free of atmospheric absorption and reflection effects (lossless data). Thus 

PSD'" = 10 logjg [?>(n)] (4b) 
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The selected relation for lossless flap noise spectra given by equations (1) and (4b) 
is based on EBF experimental evidence, acoustic noise theory, and spectral analysis 
theory. It has been shown that spectral data taken over a range of exhaust velocities for 
numerous over-the-wing and under -the- wing EBF configurations correlate when plotted 
in terms of the PSO"^ function (refs. 7 to 9). Furthermore, the functional form of equa- 
tion (4a) was used to correlate noise from jets (refs. 10 to 12). It was also used to cor- 
relate total sound power spectrum data taken over a range of exhaust velocities of vari- 
ous EBF model configurations (refs. 13 and 14). 


Theory and Experiments Related to Externally Blown Flap Noise 

The semiempirical method of representing the EBF sound spectral density was re- 
lated in reference 1 to certain theoretical and experimental work on noise source 
mechanisms believed to be related to EBF noise. Extensive work has been done in 
studying many EBF noise source mechanisms (refs. 5, 6, and 15 to 26). Some of these 
theoretical studies try to relate the turbulence spectral density to the sound spectral 
density. These theoretical studies require detailed information about the nature of the 
turbulence spectral density beyond that contained in reference 27 to make accurate pre- 
dictions. 

An important noise source for EBF systems is scrubbing noise generated in the tur- 
bulent boundary layer by the jet mixing region and convected along the wix^ or flap. 
Attempts have been made to relate the turbulence spectral density associated with scrub- 
bing to the noise produced (refs. 24 and 26). The necessary turbulence spectral density 
data used in reference 26 was obtained from references 28 to 30. 

The semiempirical method of representing the EBF sound pressure spectral density 
(ref. 1) was developed so that it could be related to a scrubbing noise source as de- 
scribed in references 24 and 26. The approach used was based on the theory of spectral 
analysis and realizable linear filters as described in references 31 to 36. Certain pro- 
cedures were adopted from modern control theory (ref. 37), The basic definition of 
one-third-octave spectra is that specified in reference 38. 


Development of Spectral Density Function 

The one -third -octave sound pressure level spectrum approximates the mean-square 
pressure spectral density of a real stationary random process. The same information 
that is present in the spectral density is present in the covariance function of a station- 
ary real random process. The equations expressing these relations are stated in appen- 
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dix B. These equations are based on spectral analysis theory. To represent a mean- 
square pressure spectrum oi a real random process these functions must be nonnegative, 
integrable, and an even function of frequency. 

The spectral density of many real random processes can be represented by a certain 
general function form. This form, called a rational function, is described in appen- 
dix C. By analogy it will be assumed that ^(f2) can be expressed as a function having 
the same basic form (see appendix C). 

A particular realization of this general form is selected. This realization involves 
the use of a function formulated from the Fourier transform of a ''damped cosine oscil- 
lation. " This function is discussed in appendix D. The resulting formula for <p{U) is 
denoted as the basic spectral density function function can be written 

as the sum of basic Strouhal response curves (appendix E). 

The basic spectral density was used in reference 1 to fit PSD'*’ data, A good fit was 
achieved with data from an under- the- wing EBF configuration and an over- the- wing con- 
figuration with attached flap flow (powered lift) and with unattached flap flow (conven- 
tional lift). Reference 1 also shows the basic spectral density function to be related to 
the theory of scrubbing noise as described in references 24 and 26. 

Although the basic spectral density function satisfies some of the requirements, a 
function that uses less computer time can be formulated from this expression. (See 
appendix F, ) The new function is empirical and since it is not an even function of fre- 
quency it is not a spectral density. However, since it has the other properties of the 
basic spectral density, the new functional form of c^(J2), denoted as the simplified spec- 
tral density function curve fit the data discussed herein. This 

function can be written as the sum of simplified Strouhal response curves (appendix F). 
Reference 1 shows that this function was also satisfactory in fitting many EBF sets of 
data. 

The simplified lossless spectral density function is 





(FI) 


where is a normalization parameter determined by the requirement that the integral 
from zero to infinity of the spectral density is unity, a^: i = 1, . . . , and bj^: k = 1, 

• * * > complex parameters with negative real parts and positive imaginary parts, 

a(n) is a polynomial equation of order and b(fi) is a pol 3 momial equation of order 
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Thus, the approach taken to represent lossless flap noise spectra is based on the 
equation obtained by substituting equation (FI) into equation (1). The resulting equation 
can be fitted to one-third- octave lossless SPL spectra with the proper selection of the 
OASPL, a^, and bj^ parameters. After developing the reflection correction, the cal- 
culation method used to select these parameters will be discussed. 


REFLECTION CORRECTION 

It has been assumed that the measured spectrum can, in part, be represented as 
the sum of a lossless spectrum and a ground reflectance. The function used to repre- 
sent the lossless spectrum was discussed in the last section. In this section the equation 
used to represent the reflectance will be derived. 

The effect of interference due to reflections is strongly dependent on test facility 
geometry. The analysis of these effects is based on the EBF test facility geometry of 
reference 2. Figure 1(a) shows the source- microphone geometry with respect to the 
ground plane. The wing flap system is shown in figure 1(b), and the microphone array 
in figure 1(c). The basic concepts used in the simple mathematical model that will be 
described for this setup are presented next to reveal some of the complications hidden 
by the approach used. 

The mathematical model deals with the reflection problem even though a good physi- 
cal description of it is lacking, be cause the data were measured in one- third- octave bands 
and thus provide only limited information representing an average over each band. Only 
one of several possible interpretations of the mathematical model will be considered in 
the analysis. The working hypothesis will be that the noise can be treated as if it were 
from a vertical (spanwise) arrangement of independent sources of equal strength at the 
trailing edge of the third flap (fig. 1(a)). The choice of this physical model to represent 
the measured effects does not indicate that this model is the best physical description of 
the true noise source. The EBF data in references 8, 13, and 14 indicate that leading- 
edge noise, trailing edge noise, flap surface noise, and redirection of jet noise are all 
possible contributors. Opinions on the dominant noise source differ. 

The effect of reflections are dealt with using an acoustic ray theory. But the ap- 
proach used differs in that the reflected ray is assumed to be related to the direct ray 
by a convolution integral of the reflecting surface's impulse response with the direct 
pressure wave. The Laplace transform of this impulse response is called the reflected- 
ray transfer function. 
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Model and Assumptions 


The interference due to reflections causes cancellations and reinforcements to occur 
in certain one- third-octave SPL bands. In the example considered, the distance from 
the source to the microphone is large compared with the source or microphone heights 
above the ground (fig, 1(c)), Thus, to a large extent, changes in the longitudinal loca- 
tion of the vertical array of sources does not shift the one-third- octave band in which a 
cancellation or reinforcement occurs. Hence, it is assumed sufficient to use, for these 
calculations, as the distance between source array and microphone the radius measured 
from the center of the jet nozzle exit to the microphone. 

To simplify the model the following assumptions are made: (1) The reflected-ray 
transfer function is the same for each of the equal independent sources in the vertical 
array. (2) At different microphone angles the reflected-ray transfer function changes. 
For a model of N independent sources of equal strength at the trailing edge of the third 
flap, this could be interpreted as being due to the reflected rays from the trailing edge 
of the third flap traveling more than one path to get to a particular microphone. The 
model having N independent sources of equal strength at the trailing edge of the third 
flap will be called the multiple reflection model. Because of the large size of the wing 
flap system (fig. 1(b)), the wing and flaps are important factors (along with the ground) 
in producing the observed reflection effects. 

The vertical displacement distances of the selected sources from the jet axis for 
each of the sources is determined by the test configuration (fig. 1), This distance yields 
the one- third- octave range over which the calculated cancellations and reinforcements 
extend. As previously mentioned, the multiple reflection model is not sensitive to 
changes in path length changes due to reflections from the wing-flap system. (It is pos- 
sible that the variations of the reflected-ray transfer function are actually due to varia- 
tions in the so\irce strength in the vertical source plane or to the interaction of the 
various possible noise sources, producing different reflection effects at different angles 
from the flap. ) 

A derivation of the reflectance that can be used to represent the effect of ground to 
wing- flap reflections on the sound produced by the wing- flap system interacting with the 
coannular jet is given next. The model assumes that the source is distributed over a 
region and approximates this by considering the rays from N independent point sources. 
This assumption was used by P. Thomas to study acoustic interference of the noise pro- 
duced by a single jet due to reflections from a plane (ref. 39). An alternative treatment 
for a distributed source using the wave equation is given in reference 40. The ground 
reflection effect on a single point source has been treated in references 41 to 44. 
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Derivation of Equations 


The noise sources, flap, and microphone geometry with respect to the groimd plane 
are shown in figure 1(a). Each of the N independent sources is assumed to produce a 
direct ray that travels a distance . All the sources are assumed to be at a fixed dis- 
tance r from all the microphones. For the n” source at a height hg(n), 

1/2 

Sj(n) = |r^ + [hg(n) - hjj (5) 

where h^ is the microphone height. Each source, as previously mentioned, is also 
hypothesized to produce multiple reflected rays that travel distances The reflec- 
tions the rays undergo vary with microphone position because different multiple reflec- 
tions are undergone by the reflected rays that reach each microphone. 

Calculations must only be accurate enough to locate cancellations and reinforce- 
ments in the proper one -third -octave band. Thus the assumption can be made that the 
path length difference between the direct -ray path and reflected -ray path from the source 
to the wing flap to the ground to the microphone is nearly equal to the path length differ - 
ence between the direct -ray path and reflected -ray path from the source to the ground 
to the microphone. This reflected -ray path distance is, for the n*^ source 

S 2 (n) = P + [hg(n) + hp} ^ (6) 


The sound reaching the microphone is assumed to consist of the sound that travels 
the direct path and the sound that travels the reflected path. Since the reflected path is 
longer than the direct path, the reflected sound reaches the microphone after a time 
delay t where 


(where c^ is the velocity of sound) relative to the direct sound that left at the same 
time. Thus if the direct wave reaching the microphone was emitted at a time t, the 
reflected wave reaching the microphone at the same time is due to one emitted at an 
earlier time t - r. 

The effect of the multiple reflections on the reflected wave is related to the direct 
wave as follows: Let p(t) be the direct normalized pressure wave incident on the micro- 


9 



phone. Then for the source the reflected wave from the wing-flap and ground sur- 
faces can be related to the direct wave, with the time delay and the distance ratio 
1/z accounted for, by the convolution integral of the reflecting surface's impulse re- 
sponse with the direct pressure wave. Thus 

— p;(t - T_) = — j g(T - T )p(t - t) dT 
Z " z„ "^0 

n n 


where 


So(n') 
z = 

^ Sj(n) 

Removing the inverse distance and the time -delay effects results in a definition of the 
effect of the reflecting surfaces: 


p;(t) = / g(r)p(t - t) dr 
n ^0 


( 8 ) 


where g(r) is the impulse of the reflecting surfaces. The reflected-ray transfer function 
G(s) is defined as the Laplace transform of the impulse response g: 


or 



g(t - 


g(r)e~^® dr 


(9) 


This complex reflected-ray transfer function can also be written in terms of its phase 5 
and amplitude Q as 

G(s) = (10) 

The observed normalized pressure at a distance r will be 
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1 


( 11 ) 


where 


N ^ 

p(t, T^) = ^ P„(t) + V p;,,(t - 

n=l Z— / 

n'=l 


r , = 5 , 
n'n n'n 


[s2(n') - Sj(nj| 


z , = 5 , 

n'n n'n 


S2(n') 

Si(n) 


( 12 ) 


(13) 


The variation of the heights of the sources, relative to the radial distance is as- 
sumed small; thus the value of will be replaced by a mean value z. Also 

will be denoted by r . 

The resultant normalized mean- square pressure spectrum density at the micro- 
phone will be found by a two-step procedure. First the autocorrelation of the signal will 
be found. Then the Laplace transform of the autocorrelation will be taken. This two- 
step procedure must be used to obtain the frequency response, because the noise is as- 
sumed to be produced by a stationary independent random process. 

The autocorrelation of the relative pressure signal received by the microphone is 

Rjjjjd) = E|p(t, T^), P(t + I, Tjj)j (14) 


where E{ } denotes taking an ensemble average (see appendix B). Substituting equa- 
tion (11) into equation (14) yields 



Since each source is assumed to be independent, the calculation can proceed usir^ one 
source, and the result can be summed to obtain the effect of N sources. For the n ^ 
individual source 


U 





JO -E^|p„(t)+^p;(t- tJ 




IL' 


nj !_' 


Pn(t-^l)-^-P;(t- Tn 4 - 0 _ 


ii 


(16) 


This equation when expanded becomes 




- ^{Pn(‘ - ""n^Piift + E{p;^(t - Tj^)p;(t - + ?)] 


= V?) J 


^{pn(t)Pn(t ■ '^n 


The Laplace transform of the autocorrelation ^(O is the spectral density of the 


source Hj^(s) with reflection effects included: 




(18) 


Substituting equation (8) and (17) into equation (18), we obtain 


Hn<®) = fj I^ot.n'^)®'^' 


-7 I g(^l)Pn(t - -^l) g(T2)P„(t + I - Tg) dTgJ d| 

z 

Z i ®‘^®E|p^(t) g(T)Pjj(t - + I - T) dTj d| 

^ I Pn<‘ + ^ + "^n)} 
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The first term on the right of equation (19) is the source spectrum of the n source 




( 20 ) 


The second term on the right hand side of equation (19) is evaluated next. Using equa- 
tion (9) and (20) 

jT g(Ti)p„(t - Tj) dTjjT g(T2)p„(t - Tg + I) dTjj d« 

J /'OO ^OO 

0 g(T2)R^(| - Tj + Tj) dTg d? 

-T,s r<=^ TgS /•» _ 

= j g(Ti)e dTj / g(T 2 )e di 2 J e R„(a) da 


= G(s)G*(s)Sj^(s) 


(21) 


The third and fourth terms on the right hand side of equation (19) are also evaluated 
usir^ equations (9) and (20). Thus 


I e'^®Efpn<‘) £ S(T)p^(t - + I - T) dr' 





J g(r) dr Rj^(^ - - r) di 


e"^®R (a) d<7 f e^®g(r) dr 
n ✓A 


= e " G* (s)S„(s) 


( 22 ) 
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and 


jT e’^®E|pj^(t + ? + T^) ^ g(T)Pj^(t - t) drj 


^<50 ^OO 

= 1 1 g(T) dT + T) 


'0 ‘'0 




-T S r®® 

= e ^ / e‘^®R (a) da / e‘^®g(T) dr 

«/q n */n 


= e ^"*G(s)Sj^(s) (23) 

Substituting equations (20) to (23) into equation (19), the spectrum due to a single 
source with reflections included is 


H„(s) = SJs) + -!■ G(s)G» (s)S„(s) + 1 .;G(s)e 

n n 2 " z 


-V 


G(s)e 


-V 


S„(s) 


(24) 


The normalized mean- square pressure spectrum density due to N independent sources 
is then 

H(s) = £ H„(s) 
n=l 


N 


N 


£ S„(s) + G(s)G* (s) S„(s) . i 

n=l n=l ^ 



<G(s)e ” + 


G(s)e ^ 


■Sn(s) 


n=l 


(25) 


The relative mean- square pressure spectrum density of the source without reflection 
effects is then 


N 


n=l 


(26) 


14 



Assume that all sources are of equal strength. Then 


Sj,(s) = NSjj(s) 


( 27 ) 


Using equation (10) 


-T^S 

G(s)e ^ 



= Q(s) 


' -T S+j6(s) 
e 


+ e 


= 2Q(s) cos [^TrfT^ - S(s)j 


(28) 


In terms of a one -third -octave spectrum, equation (28) can be easily evaluated if 
Q(s) and 6(s) are assumed to be constant in any one -third -octave band using equations 
(Bl) to (B7). Equation (28) is evaluated as foUowsr 



0^ 

27TTn 


|sin - 



sin 


[2?rf; 


a^n 



(29) 


Using the trigonometric identity 


sin 




- sin co„T„ =2 cos 


H 


a n 


Tj^sin 


n 


(30) 


and equations (B6) and (B7), we obtain from equation (29) 



Q(s) cos - 6(s)J df = Af^Q(f^) cos 




sin (j7 Af^Tj^) 




(31) 


If the source spectrum Sp,(s) is assumed to be constant over any one-third-octave 
band the measured one -third -octave spectrum (neglecting atmospheric attenuation) can 
be obtained by substituting equation (25) into equation (B2) and using equations (10), (27), 
and (31). Thus, 
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r 


Q(fc) 




L 



r -I sin {n Af t ) 

cos [2rf^T^ - e{f^)] 




n=l 


J 

(32) 


where the notation has been changed in that since s = parameters hitherto repre- 

sented as a function of s will henceforth be represented as functions of f 

This equation represents the output one-third-octave spectrum H(f^) Af^ as the 
product of the input one -third- octave spectrum Sp(f^) Af^ and a reflectance factor 
T(fg)p, where 


|T(f^)|^ = 1 +i-Q(f 

'Cl ^2 ^ N z 



n=l 


P sin (JT Af T ) 

cos [2rfj.T„ - 5(f^)] 


(33) 




Thus 

H(f,) Af^ = Sp(f^)|T(f^)f Af^ (34) 

Equation (34) can be interpreted in terms of well-defined single input, which is from a 
stationary random process having a one-third- octave spectrum Sp(f^) Af^ and an output 
that belongs to a stationary random process having a one-third- octave spectrum 
H(f,) A(^. 


Reflected-Ray Transfer Function 

j6(fp) 

The reflected-ray transfer function G(f^) (where G(f^) = Q(f^)e , eq. (10)) 

must be specified to calculate the reflectance |T(f ) (eq. (33)), If the reflected-ray 
transfer function depended only on the acoustic impedance of the ground surface, its 
value might be measured directly, using the method discussed in reference 45, or calcu- 
lated based on the assumption that the ground acoustical properties are similar to the 
measured properties of some material (refs. 46 to 49). But the problem of evaluating 
the reflecting properties of a ground surface and using the results with an extended 
source such as is produced by an engine or nozzle is so complicated that common exper- 
imental solutions to the problem are either to apply a sound absorbing material to the 
reflecting surface (refs. 50 and 51) or to use flush mounted microphones (refs. 52 
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and 53). The problem becomes complicated when the source is an externally blown flap, 
because, for short wavelengths reflections from the wing -flap and ground, surfaces may 
become randomized such that the interference effect becomes to some extent indeter- 
minate. Then the reflected-ray transfer function may be determined by a scattering of 
the incident wave . 

Specific factors important in evaluating the reflected-ray transfer function are 
(1) the variation of the surface impedance with moisture in the soil and vegetable roots 
in the soil (ref. 45) and (2) the variation of the reflection coefficient Cj^ with the angle 
of incidence The reflection coefficient is related to the surface impedance and 

'’inc ’’y 


CR(f. ^inc) = 


Z(f) cos 9.^^ - pCp 
Z(f) cos + pc^ 


(35) 


where Z(f) is the surface impedance (in general, a complex quantity dependent on fre- 
quency) in acoustical ohms times square meter (also called rayl), p is the density of 
air, c^ is the velocity of sound in air (pc^ = 416 J2*m^ (MKS rayles)), and is the 

angle between the incident ray and the normal to the surface (ref. 44). The absorption 
coefficient of the ground is related to the reflection coefficient Qf(0) by 

|Cp(f,0)p = 1 - a(f,e) (36) 

The reflected-ray transfer function G(f^) would be equal to the reflection coefficient 
Co if the reflected-ray transfer function depended only on the ground impedance. Thus, 

rC 

G(fc) =Cjj(fc,fl) (37) 

and the reflected-ray transfer function G(f ) is related to the absorption coefficient 

|G(fc)r =[Q«c^ = 1 (38) 

The reflected-ray transfer function G(f ) is assumed to be independent of frequency. 
This assumption is made, first, because impedance tube measurements of the normal 
absorption coefficient = a(f,90*^) of sharp sand over a frequency range from 200 to 
1000 hertz for a moisture content ranging from 0 to 10 percent indicated that the normal 
absorption coefficient was not strongly dependent on frequency at a given moisture 
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content (ref. 45), However, the normal absorption coefficient changed with moisture 
content. Thus equation (38) indicates that under some conditions certain materials have 
reflected-ray transfer function magnitudes that do not change with frequency but do 
change with m.oisture content. The reflected-ray transfer function is assumed to be in- 
dependent of frequency, secondly, because of a lack of ground surface impedance data 
at the site of theEBF turbofan ermine experiment and because of the possibility of scat- 
tering being an important factor in the multiple reflection model introduced to deal with 
the problems created by an extended EBF source. Thus 

G(fj) = Qei® (39) 

where the constants Q and 6 are now frequency independent. 

CALCULATED SOUND PRESSURE LEVEL 

The calculated sound pressure level equation used herein is formed by first com- 
bining the one -third-octave lossless sound- pressure equation (eq. (1)) with a spectral 
density equation (eq. (FI)) and equations for the reflection correction (eqs, (33) and 
(39)). Thus 

SPL’’ai&c»)] = 10 

= SPL+(0+ 10 logic |T[fc(2)]|^ 
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where OASPL^ is a parameter representing the lossless overall sound pressure level, 
is spectral density normalization parameter, a.:i = 1, . . and bj^rk = 1, 

. . • > Njj are complex parameters that aid in determining the shape of the spectral 
density curve, Q is a frequency independent parameter that determines the amount the 
reflected wave interferes with direct wave, and 6 is a frequency independent parameter 
indicating the phase angle change in the reflected wave relative to the direct wave. This 
spectrum must then be corrected for atmospheric attenuation. The correction used in 
this report is based on that in reference 54. Let A^Q^(f^,r) represent this correction; 
then 




- SPL'^(Z) + 10 logjQ 



+ A 


tot 


fcU),r 


(41) 


Equation (41) will be used to represent the measured spectrum. The calculation 
method used to identify parameters that will enable equation (41) to fit a one- third- 
octave band SPL spectrum will be discussed next. 


CALCULATION METHOD 

A computer method is used (1) to find an empirical curve that represents the meas- 
ured SPL corrected for reflection effects and atmospheric absorption and (2) to conduct 
a numerical evaluation of the effect of reflections on the measured SPL. This is 
achieved by using the computer to determine parameters such that a calculated reflec- 
tion correction and a calculated lossless SPL spectrum are adjusted to give a minimum 
least- square error when compared with a measured SPL spectrum after taking into 
account atmospheric absorption. 

The computer program is initiated by first selecting parameters that define a math- 
ematical model for the ground reflection effect and for the sound spectrum and then, the 
initial values of the parameters, which the computer adjusts. The adjusted parameters 
are subject to certain constraints. These constraints are used to define cost and penalty 
functions (refs. 55 and 56) to assure that the final parameters satisfy the constraints. 


Cost Function 

The computer program is used to calculate the SPL^^j parameters (a., OASPL^, 
Q, 6). The program is designed to minimize a cost function by adjusting the 
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parameters. Actually, the computer program adjusts the parameters to minimize a 
quantity which is the sum of the cost function and a penalty function 

The cost function is the weighted sum of the squares of the differences between 
the measured SPL_ and calculated SPL^^i (eqs. (40) and (41)) sound pressure level 

HI 

at each frequency fp(0* Thus, 


. c 

Ce = Z “e{sPL„[fc(0] - SPL,^l[yZ)]| 


(42) 


The total cost function is thus defined as 


2 

Ctot ‘ ®"’^cai[fc<^)]} - Sen 


(43) 


The penalty function, Cp^^, is discussed next. 

As mentioned previously the parameters are subject to constraints. (See eq. (2) 
for the constraint on OASPLp and eq. (5) for constraints on a. and bj^. ) To insure 
compliance with these constraints a ’’penalty function” Cp^^^ is added to the cost func- 
tion. The penalty function increases the cost when an associated constraint is violated. 
Since the computer iteratively adjusts the parameters to minimize the total cost function, 
parameters violating a constraint are not reached by the adjustment process. 

The constraints arise as follows. The parameters in equation (5) have limits. They 
must be chosen so that the poles bj^ and the zeroes a^ of the function <p (Q.) have nega- 
tive real parts and positive imaginary parts. Also, the values of the parameter 
OASPL^ must satisfy equation (2); thus 


OASPL+ = 10 log^Q 



(44) 


The penalty function has the following characteristics. If the real part of a root of 
a(fl) or b(fl) is less than zero, or the imaginary part of a root of a(fl) or b(fl) is 
greater than zero, the individual penalty is zero. If the real part of a root of a(fl) or 
b(fl) is greater than or equal to zero or the imaginary part of a root of a(l^) or b(J2) is 
less than or equal to zero, the individual penalty is equal to the square of the real part 
of the root or the imaginary part of the root. The total penalty is the sum of the indi- 
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vidual penalties for all the roots plus the weighted sum of the squares of the differences 
of the parameter value of the reflection corrected unattenuated overall sound pressure 
level OASPLp and the calculated value of the reflection corrected unattenuated overall 
sound pressure level OASPL^^j (based on eq. (43)). Thus, 

“a 

C = 6 

pen / 

i=:l k=l i=l 

k=l 

where 



+ iS^OASPL^^ - OASPL+ j (45) 


N. 


M. 


^(i)|'/?ea.j^ + J ^ 6^(i)j:j?/» a^j‘ 



a. S: 0 
Otherwise 

Otherwise 
a. < 0 
Otherwise 
jW bj^ < 0 
Otherwise 


(46) 


(47) 


(48) 


(49) 


Optimization 

The problem of minimizing a nonlinear function like the cost function given by equa- 
tion (43) can be solved by search techniques. These techniques consist of systematic 
procedures for varying the cost function parameters until a minimum value of the cost 
function is found. The basic methods available are discussed in references 55 and 56. 
To provide flexibility in the choice of the order of equation (15) (i.e. , the number of 
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terms and Nj^), and to accommodate a penalty function, a search technique that 
does not require evaluation of derivatives was chosen. The search technique was that of 
Powell (ref. 57). The computer program was adapted from reference 58. The optimi- 
zation procedure has previously been applied successfully to a similar identification 
problem (ref. 59). 

For modeling acoustic data choosing the number of terms in the linear system 
power density spectrum function (i.e. , the model order) is a serious problem. The 
selection of a low- order model limits the frequency range over which the function re- 
sembles the data. The selection of a high-order model causes the function to adjust to 
all the variations in the measirred data including those due to the reflection effects. The 
effect of adjusting the model order (i.e., changing the number of poles or the num- 
ber of zeroes M„ in eq. (40)) on the minimum residual error provides some guide as 
to the appropriateness of the model. Thus the variation of the minimum cost obtained 
for different values of and helps determine the best values to use. Various 
schemes besides the test of the residual error have been proposed to help select the 
model order for identification problems (ref. 60). 

The number of poles selected for the model was one. The number of zeroes 

was three. These values were selected specifically for the data discussed in the 
application section of this report. 


RESULTS 

APPLICATION TO MEASURED DATA 

The method for finding the function and the parameters Q and 6 was 

applied to data from a large scale externally blown flap. The scope of this report does 
not include interpretation of these data, since they are interpreted and discussed in 
reference 1. 

An N independent point source model consisting of five point sources of equal 
strength was selected to represent the data which the procedure was applied to fit. One 
point source was chosen to be on the engine axis at the nozzle exit. Two other point 
sources were assumed to be 36. 6 centimeters (1.2 ft) vertically above and below this 
point, and two more point sources were assumed to be vertically displaced 73.2 centi- 
meters (2, 4 ft) above and below the engine centerline. Calculations show that the use 
of more than five point sources does not appreciably affect the magnitude of the cancel- 
lations and reinforcements so long as the total maximum vertical displacement is the 
same. 
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Test Conditions 


The test data to be evaluated were obtained using a coannular nozzle. The flaps 
were set at 0°-2 0^-40° flap angle pattern, meaning that the leading flap was not de- 
flected, the trailing edge of second flap was deflected 20° from the wing reference chord 
line, and the trailing edge of the third flap was deflected 40° This is referred to as the 
takeoff flap condition. A schematic diagram of the flap configuration is shown in fig- 
ure 1(b). The nozzle axis was at a 5° angle with respect to the reference chord line of 
the wing. More characteristics of the test configuration are given in table 1. 

The microphones (fig. 1(c)) were placed at a radius of 30. 48 meters (100 ft) and 
were at the engine centerline level, which was 2. 74 meters (9 ft) above the ground. 

This is only the nominal height, since the ground at the site sloped slightly. The micro- 
phone angles from the engine inlet considered in this report are 0°, 20°, 40°, 60°, 70°, 
80°, 90°, 100°, 120°, and 140°. Most emphasis is placed on the 100° microphone loca- 
tion since data from three microphones at different heights (1.37, 2. 74, and 5. 49 m 
(4. 5, 9, and 18 ft)) were measured at this position. More details of the test and turbo- 
fan engine are available in references 1 to 3. 


Test Data 

Representative one -third- octave SPL plots of the data from this externally blown 
flap (EBF) turbofan engine configuration are shown in figure 2. The SPL data selected 
are from microphones at 20°, 90°, 100°, and 120° for the maximum velocity test. 

The plots show that at frequencies above 4000 hertz engine machinery noise sources 
are probably having an effect on the data. This assumption is also based on other EBF 
data presented in references 13 and 14. Therefore, only the SPL data at frequencies 
from 50 to 3150 hertz are used in the following discussion. 

From equation (40) the first cancellation due to the direct and reflected waves being 
180° out of phase for a single point source and a hard surface (Q = 1, ^ = 0), would occur 
at a frequency 


f. = (50) 

2(S2 - Sj) 

and the second reinforcement due to the direct and reflected wave being in phase would 
occur at 2f^. For the test facility shown in figure 1, f^. is within the 3 15- hertz one- 
third- octave band when the microphones are at 2. 74 meters (Oft). The first cancellation 
is indicated on the SPL plots by and the second reinforcement by R2 (R^ is at 
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0 hertz). For microphones at 1. 37 meters (4. 5 ft) ^he first cancellation is within the 
630 -hertz band, and the second reinforcement is within the 12 50 -hertz band. For mi- 
crophones at 5. 48 meters (18 ft) the first cancellation is within the 160 -hertz band, and 
the second reinforcement is within the 315-hertz band. 

Examination of the SPL at the 20° microphone shown in figure 2(a) indicates a can- 
cellation at 315 hertz. A cancellation at 315 hertz islalso apparent in the SPL from the 
90° microphone shown in figure 2(b). No cancellations are indicated in the SPL at the 
100° microphone shown in figure 2(c) or the 120° microphone shown in figure 2(d). If 
the ground is assumed uniform, the ground impedance! concept would give equal interfer- 
ence effects at all angles. The data plots indicate that this is not the case. The concept 
that a reflected- ray transfer function could vary with T^iew angle and cause the different 
interference effects at different microphones is an explanation for these measurements. 
The multiple reflection model is emphasized herein, sijnce it is known that the relative 
position of the flaps to the hypothesized N independent^^ sources on the trailing edge of 
the third flap changes with the angle of observation, while the noise source azimuthal 
directivity variation, the nature of the dominate noise source and the interaction of the 
various noise sources are not known currently. | 

Effective diameter and effective velocity are used in the conversion of SPL'*’ to the 
PSD'*’. The effective diameter is calculated as follows: 




^ ^^fan ^core^ 


1/2 


The effective velocity is calculated from 




/A V® +A V® 
^'^fan ^fan ^ '^core^ core 

■^fan ^core 


(51) 


(52) 


where the velocities used are obtained by isentropic expansions from total to ambient 
conditions. 

These relations (eqs. (51) and (52)) are based on the data described in references 7 
and 8. For the data of reference 8, these relations applied reasonably well to the 
simple nozzle -flap configuration tested using a convergent and a coannular nozzle. The 
tested configuration of reference 8 was based on one of the double slotted external flow 
jet flap configurations developed by the NASA Langley Research Center. Also in refer- 
ence 8, a scaling procedure is described that can be used to estimate full-scale blown 
flap noise from small scale data. 
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Calculation Procedure 


Figure 3 shows a diagram of the computer calculation method. The input informa- 
tion consists of the following: (1) the measured sound spectrum SPL^, (2) the con- 
stants defining the spectrum density equation order (M^ and N^j), (3) the constants used 
to convert from Strouhal number to frequency (D^ and Vg), (4) the constants defining 
the reflection model (N, z, r^:n = 1, . . . , N), (5) the initial parameter values 
(OASPLp; a.:i 0, . . . , M^; = 1 , , N^; Q; 5), and (6) the atmospheric atten- 

uation correction. The computer program calculates at the start of each iteration step 
the following quantities using the current set of parameter values: (1) the reflectance 
(|T(fg)|^), (2) the mean- square-pressure spectrum density level (3) the SPL 

spectrum without reflection and atmospheric attenuation effects SPL^^^p (4) the SPL 
spectrum with atmospheric attenuation and reflection effects SPL^^^P Next, the cost 

function C is calculated from equation (42). Then the penalty function C is calcu- 
^ P 

lated using equations (43) to (48). Last, these two quantities are summed to obtain the 

total cost Ctof The new total cost is compared with previous values of the total cost 
to see if another iteration is necessary. The calculations may stop at this point or con- 
tinue. If the process is to continue, the new total cost and previous values of the total 
cost are used by the computer to adjust the current parameter values and to obtain new 
parameter values that will reduce the total cost. These new values are used to initiate 
a new iteration. 

If another iteration is not needed, the procedure terminates. The current values of 
the parameters I T(fg)|^, P^D^^p SPL^^p and SPL^^^ are then available. Besides 
these quantities, a mean square pressure spectrum density level based on the measured 
SPLjj^ and the calculated reflectance and atmospheric attenuation PSD^ is also cal- 
culated. This can be compared with to check the accuracy of the curve fit and 

the ability of the reflectance to account for interference effects. 

The process is thus seen to be based on comparing a measured sound spectrum 
SPLj^ with a hypothetical calculated sound spectrum SPL^^P The hypothetical calcu- 
lated sound spectrum SPL^^^j^ is based on two calculated curves: (1) a cvirve 

and (2) a |T(fg)|^ curve. These curves are assumed to be independent. The final curve 
shapes are achieved by the iteration procedure are assumed to model the actual physical 
values of these quantities. 

As mentioned previously one of the first steps is to select the PSD^^^ model and 
the values of the PSO^^I parameters. This is covered in detail in appendix F. The 
reflection model parameters and z are based on the test site geometry. Numerical 
experimentation indicates that increasing N beyond five does not change the reflectance 
if the distance between the furthest source and the nozzle exit centerline position is fixed 
for fixed z, Q, and 5. The value of Q is between zero and one and the initial 6 was 
selected to be 2ir. 
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Comparisons 


Initial and final iteration values of SPL^^j, SPL^g^p |T(f^)| , and are 

shown in figure 4. Figure 4(a) shows the initial and final hypothetical SPL^^^ curves. 
Also for comparison, data for the input sound spectrum SPL^ are shown. Although 
the initial hypothetical spectrum does not agree with measured spectrum, the computer 
program was able to adjust the parameters so that the final calculated spectrum is in 
good agreement. The largest error occurs near 800 hertz. 

Figure 4(b) shows the initial and final SPL^^^ curves. The final SPL^g^^ curve 
has a peak level of 107. 16 decibels at 80 hertz. 

The initial and final reflectance level curves are compared in figure 4(c); the final 
value of Q was 0. 48. Figure 4(d) shows that the final curve had a peak level 

value of 4. 22 decibels at a Strouhal number of 0. 244. At high Strouhal numbers both the 
initial and final curves have similar slopes. 

Figure 4(e) shows PSD"^ based on correcting the measured data for atmospheric 
attenuation effects and for reflectance effects. The reflectance effects are corrected 
using the hypothetical reflectance level curve. The resulting PSD^ curve is compared 
with the final PSD"^ curve based on the curve fit PSD^gp Again, the error is largest 
near a Strouhal number of 31, which corresponds to 800 hertz. The final calculated 
PSD^ai curve is in good agreement with the curve based on measured data. 

In appendix G characteristics of the rational function PSD”^ and trends of the curve 
fit values of the OASPL^ parameter are discussed. 


REFLECTION CORRECTION 

Two aspects of the reflection correction will be covered. The first deals with 

ability of the method to correct or decrease the cancellations and reinforcements in the 

SPL and the second with the reflected ray transfer function, 
m 


Application of Reflectance 

The concept to be illustrated is the existence of a hypothetical PSD'*’ curve and a 
hypothetical reflected ray -transfer impedance at a given angle and/or a given test con- 
dition that can be used to calculate the SPL that would be measured by microphones at 
different heights and distances from the source. The assumption is made that directivity 
and angle of incidence effects are secondary and may be neglected. The results obtained 
using the procedure with data measured at the 100° microphone at the highest exit veloc- 
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ity test condition (V^ = 242.3 m/sec (795 ft/sec)) are shown in figure 5 for three micro- 
phone heights. Figure 5 shows a plot of the calculated unaffected sound pressure level 
compared with the sound pressure level including reflection effects and atmos- 
pheric attenuation and the measured sound pressure level SPL^. 

The agreement between the measured data and the calculated SPL^^^j is good. The 
cancellations and reinforcements occur at nearly the same frequency and are of similar 
size in both measured and calculated SPL curves. The error is of the order of 1 deci- 
bel at the frequency of maximum cancellation. Moreover, the calculated SPL"** should 
provide a satisfactory estimate of the actual SPL^ since the SPL^^ agrees with the 
SPL^^j at different heights. 


Reflected- Ray Transfer Function 


For the previous case the value of the reflected- ray transfer function coefficients 
after the cost function was minimized were Q = 0. 1 and 6 =5.94 radians (+340°). To 
obtain some idea of how a change in the magnitude of the reflected- ray transfer function 
Q affects the data, the case shown in figure 5 can be compared with the one shown in 
figure 6. Figure 6 presents the data and results for a case having a Q value of 0. 5. 
This is the largest value of Q evaluated by the procedure in the study of the data. The 
maximum cancellation is shown to be about 5 decibels in figure 6 for the Q equal to 
0. 5. This can be compared with a maximum cancellation of about 2 decibels for the 
data shown in figure 5, in which Q = 0. 1. 

Values of Q and 6 are presented in table n. Since the value of Q indicates the 
magnitude of the effect of the reflected- ray transfer function on the cancellations and 
reinforcements, it will be discussed in some detail. Figure 7 shows a plot of Q with 
angle over the range of velocities used in this study. The general trend shown is that Q 
is largest from 20^ to 60^ and smallest from 100^ to 120^. 

To explain this, consider that in previous treatments of reflections (refs. 41 to 44), 
the reflected -ray transfer impedance could be considered to be caused by the ground 
surface impedance, since only the reflection from the ground plane occurred. The re- 
flection coefficient calculated in reference 44 is equivalent to the reflected-ray 
transfer function G when only ground reflection is treated. Thus if the ground imped- 
ance is assumed to be independent of frequency of equation (35), 


Z cos 0.^^ - pcg ' 


R 


(53) 
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The ground material at the test site at which the EBF data were obtained is uniform in 
all directions about the wing flap system, and thus the ground impedance at any angle 
should yield the same value for the reflection coefficient and produce the same reflection 
effect at each microphone. Since the reflection coefficient depends only on the surface, 
no change in the effects it produces should occur with changes in microphone angle. 

The observed variation of Q thus seems unrelated to the usual reflection coefficient. 
One hypothesis that explains the observed variation is that multiple reflections occur due 
to the presence of the flaps. If this is assumed, the multiple reflections that occur to a 
ray as it travels from the source to a microphone would be different at each microphone. 
But, due to the geometry they would be similar for all sources in the same vertical line. 
Thus reflection effects would be expected to be different at different angles. 


SUMMARY OF RESULTS 


It has been shown that a computer procedure that corrects for cancellations and re- 
inforcements and that provides for a satisfactory curve fit to sound pressure level data 
can be designed and implemented. The most important results obtained are that 

1. A rational polynomial with a first-order numerator and a third -order denomi- 
nator is a satisfactory analytic expression for the Strouhal normalized pressure spectral 
density function. 

2. The N independent source model of Thomas is satisfactory to use for correct- 
ing externally blown flap data obscured by the presence of interference effects of reflec- 
tions if it is extended to include a reflected-ray transfer function. 

The method discussed in this report has been applied successfully to engine- over- 
the-wing and engine -under -the -wing externally blown flap data (ref. 1) and is applicable 
to nozzle only data. 

A natural extension of this work would be a study of the way the parameters that de- 
termine the PSD'*’ (lossless Strouhal normalized mean square pressure spectral density 
level) function vary with angle and velocity. A set of data that provides more informa- 
tion on the low-frequency roll-off characteristics of the PSD"^ should be used. The 
PSD'*’ function parameters found in this report were not studied to determine trends with 
velocity and angle since most of the data concerned only the high-frequency roll-off por- 
tion of the PSD'*' curve. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 19, 1974, 

501-24. 
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APPENDIX A 


SYMBOLS 


A 

core 

^fan 

A(s) 

a($2) 


B(s) 

h{Q) 


^0 

^pen 

CR(f^^inc) 

^tot 

Cl, C2, Cg 


c 


E 

Ea 

E{} 

F{n,n,c) 


core area, cm^ (in. 

fan area, cm^ (in. 

* 

zero of (p (s), sec'^ 
polynomial equation for zeroes of <;^)(s) 
atmospheric attenuation, dB 
zero of <p{^), dimensionless 
polynomial equation for zeroes of (p{^) 
a complex number, normalization constant 
pole of <p{s), sec'^ 
polynomial equation for poles of cp{s) 
zero of (p{^), dimensionless 
polynomial equation for poles of 

cost function, sum of square of error between measured and calculated 
sound spectra 

constant 

penalty function 

reflection coefficient 

total cost function 

cancellation frequencies 

a complex number 

velocity of sound, m/sec 

effective diameter, m 

see eq. (Dll) 

see eq. (D12) 

see eq. (D13) 

expected value of { } 

Strouhal response function 



G(s) 

g(T) 

H(f) 


H(fc) 


hg(n) 

I 

N 

Nu 


OASPL 

OASPL 

P(t, V 

P(co) 

PSD 

PSO"*" 


P 


PSD 


cal 


pso;)' 

m 


PSD 

p(t) 


real 




frequency, Hz 

nominal lower one-third-octave band-edge frequency, Hz 

nominal upper one -third -octave band- edge frequency, Hz 

one -third- octave center frequency, Hz 

critical frequency, Hz 

frequency of first cancellation, Hz 

reflected- ray transfer function 

reflecting surface’s impulse response function 

relative mean-square pressure spectrum density per unit frequency of sig- 
nal referenced to 2x10”^ Pa 

one -third- octave relative mean -square pressure spectrum density of signal 
relative 2x10 ^ Pa 

microphone height, m (ft) 
height of n source, m 
one-third-octave band number 

number of numerator a(f2, a^^) terms of number of numerator 

{Q - terms of 

number of independent noise sources 

number of denominator bj^) terms of number of denomi- 

nator - bj^) terms of 

overall sound pressure level relative to p^, dB 
lossless OASPL parametric value, dB 
observed relative pressure, dimensionless 
mean-square pressure spectrum density 

Strouhal normalized mean-square pressure spectral density level, dB 
PSD without atmospheric absorption and ground reflection effects, dB 
PSD"^ calculated from formula, dB 
PSD''’ calculated from measured data, dB 

PSD"'’ calculated from a real random process spectrum density, dB 

-5 

relative pressure perturbation signal referenced to 2x10 Pa, dimension- 
less 
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^tot 

p'(t) 

Po 

Pr 

Q 

Rp(r) 

^ 2 * ^ 3 ’ ^4 
r 

S^{i) 

SPL 

SPL'*' 


SPLcai 

St 

S 

®2 

|T(gl^ 

Ta 

t 


u 

V 

core 



X 


total pressure, Pa 

reflected relative pressure perturbation signal referenced to p^. 
average atmospheric pressure, Pa 
reference pressure 2x10''^ Pa (2x1 0“^ M bar) 
magnitude of reflected- ray transfer function 

relative mean square pressure autocorrelation function referenced to 
2x10 Pa, dimensionless 

relative mean square pressure autocorrelation function of signal referenced 
to 2x10 Pa, dimensionless 

reinforcement frequency 

distance from source to microphone, m (ft) 

relative mean square pressure spectrum density per unit hertz of all 
sources referenced to 2x10 Pa 

one-third-octave sound pressure level reference to p^,, dB 

one -third- octave sound pressure level referenced to that would be 
measured without atmospheric attenuation and in the absence of reflect- 
ing surface, dB 

SPL calculated from an equation, dB 
SPL measured, dB 
Strouhal number, f(Dg/Vg) 
j2irf, jw 

distance from source to microphone along direct ray path, m 
distance from source to microphone along reflected ray path, m 
reflectance 

averaging time interval, sec 
time, sec 

frequency or Strouhal number ratio, dimensionless 

core velocity, m/sec 

effective velocity, m/sec (ft/sec) 

fan velocity, m/sec 

horizontal position of wing relative to engine (see fig. 1(b)), m (ft) 
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y 

z 

z 


a(f,ei„c) 

01 


I 


a 

/3 


nor 


y 

Ap 

5^(i),5i3(k) 

5c(k),6^(k) 

nn 

6(s) 

9 

^inc 

n 

P 

c) 


'^n 

^2 


vertical position of wing relative to engine (see fig. 1(b)), m (ft) 
ground surface impedance, acoustical ohm • m (mks rayl) 
mean ratio of path length reflected ray traveled to path length direct ray 
traveled for all N sources 

absorption coefficient 

weighting factor, dimensionless 

normal absorption coefficient, a(f , 90^) 

dampir^ constant 

weighting factor, dimensionless 

constant 

constant 

constant, 1/sec 

bandwidth at center frequency (0.2316 f^.), Hz 

pressure perturbation, Pa 

see eqs, (46) and (47), dimensionless 

see eqs. (48) and (49), dimensionless 

kronecker delta (n / = 0; n = = 1) 

reflected- ray transfer impedance phase factor, radians (deg) 

frequency response excitation number 

angle between mean flow direction and observation direction, deg 

acoustic ray angle of incidence, deg 

normalized time dealy (m = 

2 

density of air, l?^/m 

p 

acoustical impedance of air, 416 acoustical ohm vcT (416 mks rayl) 

function of Strouhal number such that reciprocal of Fourier transform of 

2 

damped oscillation is proportional to a(J2, c) 
time delays, sec 

acoustic delay time from source to microphone, sec 
spectral density function, dimensionless 
relative angular velocity, j2;rf = j27rSt 
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angular velocity, 2?rf, rad/sec 


(jt) 

Subscripts: 
cal calculated 

e effective 

i index of roots of numerator of rational function cp 

inc incidence 

k index of roots of denominator of rational function cp 

m measured 

n source index 

o constant 

p parametric 

pen penalty 

real real 

tot total 

Superscripts: 

* complex conjugate 

+ spectrum with reflection and atmospheric attenuation effects removed 

*’ spectrum with atmospheric attenuation effects removed 

reflected ray 
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APPENDIX B 


RELATION BETWEEN THE ONE -THIRD- OCTAVE SOUND PRESSURE LEVEL, 

SPL, AND THE PRESSURE COVARIANCE FUNCTION, Rp(T) 

The one- third-octave sound pressure level spectrum SPL(i) is related to the mean- 
square -pres sure spectral density H(w) of a stationary random process by 


where 


and 


SPL(Z) = lOlog^Q AfpU)} 

(Bl) 

H[fc(Z)] Af^(Z) = H(a>) do) 

*a 

(B2) 

f (Z) = 

(B3) 

f = 2'^/^f ^ 0.8909 f 
a c c 

(B4) 

% = ^ 1. 1225 fj. 

(B5) 

f - ff f 1^/2 ^ ^ 3- 

(B6) 

Af^ = 0.2316 f^ 

(B7) 


The mean-square pressure spectral density H(ce) will be assumed to be defined 
over the interval from zero to infinity and referenced to a pressure of 2x10 ^ pascal 
(2x10"'^ jubar). Thus it can be related to the mean-square pressure spectrum density 
defined over the interval from minus infinity to infinity as follows 


H/a,) = 


2 


(B8) 



where 0 < w < « . The mean -square pressure spectrum density P(cu) contains the 
same information as the pressure covariance function, Rp(r). The pressure covariance 
function is calculated by 

T 

iF Jt - Po] [Ptot<‘ + - Po] 

•^a “^-^a * a 



/ T 

^ Ap(t) Ap(t + t) dt 

-r 
a a 


For a stationary real random ergodic process, the time autocorrelation function of a 
sample function of a random process equals the ensemble average (statistical autocorre- 
lation function) of that process. Thus 


Rp(r) = E - p„)2} = E {(Ap)2| (B9) 

The expression E{ } denotes taking an ensemble average. However, the statistical 
autocorrelation can also be used to define the autocorrelation of a complex random pres- 
sure processes. Let Apj and Ap 2 be the random variables that refer to the possible 
values which can be assumed at the time instants t^ and t 2 , respectively, by the 
sample functions Ap(t) of a given complex random pressure process. Then the autocor- 
relation function is defined as 


Rp(tj,t2) = E(Apj Ap|) (BIO) 

where the asterisk denotes the complex conjugate. 

If the given random pressure process is stationary, the joint probability distribution 
for Apj and Ap 2 depends only on the time difference r =t^ - t 2 - Thus, 

Rp(T) = Rp(t,t - T) = E(Ap^ Ap*_^) 

Also, the indicated averages are invariant under a translation of the time origin. Now 

E(Apj Apj*,^) = E(Apj,_^^ Ap* ) = E» (Ap^, Ap*^^ (Bll) 

therefore , 

(B12) 


RpW = r;(-t) 
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Thus the autocorrelation function of a stationary real random process is an even function 
of T (ref. 32). 

The mean-square-pressure spectrum density P(w) is the Fourier transform of the 
pressure covariance function. 


-oo 

Rp(T)e‘J“’’dT 


(B13) 


where -« < a? < <» . A function, in order to represent the mean-square pressure spec- 
trum density P(a?) must be nonnegative and integrable (finite power) (ref. 32). The 
pressure spectrum density of a real random process must be an even function of fre- 
quency. This follows from equation (B13), since for a real process 

Rq(t) = R(-t) (B14) 

(ref. 32, p. 106). 
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APPENDIX C 


ORIGIN OF RATIONAL FUNCTION FORM OF (fi) 

The spectrum density of a stationary real random physically realizable process can 
in many cases be represented as a rational function of frequency f: 


\JT 2 

cp is) = = A£sl (C; 

% B(s) 

IT 

k=l 

where s is a function of frequency defined by s = j 2 ; 7 -f (where j = -y^), the zeroes of 
B(s) (i.e., the Np poles) are denoted Bj^), the zeroes of A(s) (i.e., the zeroes) 
are denoted Aj^, and A^ is a normalization constant defined such that 


/ oo 

^(f) 


df = 1 


(refs. 32 (p. 233), 33 (p. 378), 34 (p. 410), and 35 (p. 187)). 

This functional form of the spectral density can be related to a number of physical 
phenomena. One example is the spectral density of the output from a lumped- parameter 
linear electric circuit, which is excited with white noise (ref. 32, p. 105). It resembles 
measured one-dimensional isotropic turbulence energy spectra (ref. 27). This func- 
tional form can also be used to describe the spectral density of other stationary random 
processes (ref. 36), This form could also be used to represent the power spectral 
density of broad- band aerodynamic noise with and without the presence of solid boundar- 
ies (refs. 18, 21, and 23), According to reference 34, any spectra can be approximated 
sufficiently closely by a rational function of s. Since <;^)(s) is a spectral density, it has 
particular properties that imply certain restrictions on the number and location of its 
poles and zeroes. The parameters A- and must be chosen so that ^(s) is a real, 
even, nonnegative, integrable function of the frequency f whose inverse Fourier trans- 


form is a correlation function. The zeroes of A(s)(Ap and the zeroes of B(s)(Bj^) are 

constrained to have negative real parts. Also, all A. and Bj^ with nonzero real parts 

must occur in conjugate pairs. Another condition is that the degree of the numerator 

Mg must be less than the degree of the denominator N (N < N ; ref. 32, p. 233). 

P 2 p 
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The function of Strouhal number <p (fi) will be assumed to be defined by the Fourier 
transform of a function of a normalized time delay (where m = denoted as 

R(m). Thus 

P(n)= f R{M)e‘”'^dM 


where 


^ = s — = j2jrf — 
V. V, 


= ]2j7St 


In the Strouhal number domain the function of Strouhal number will be 

assumed to be represented by the following rational function of S2: 


<p{a) = 


AoTT(« - Ai) 

1 = 1 

TT - B^) 

k=l 


A(S2) 

B(f2) 


2 


where the zeroes of A(^^)(A^ and the zeroes of are constrained appropriately 

as mentioned previously in discussion of the function ^(s) defined in the frequency do- 
main and Aq is again a normalization constant. 
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APPENDIX D 


SPECTRUM OF A ^^DAMPED COSINE OSgLLATlQN* * 

Let the pressure covariance function Rp(T) be represented by a ’damped cosine 
oscillation' ' 


Then from equation (B14) the 


Rp(r)=re cos 

mean square pressure spectrum of Rp(r) is 
P(co) 

|cr(a),c)p 


(Dl) 


(D2) 


where 


o(w, c) 



c 


-“o + j^o 


(D3) 

(D4) 


and the. asterisk (* ) indicates the conjugate relation. 

The behavior of the function of co, P(w), for different values of c is more easily 
discussed if the function is written in terms of dimensionless variables. Let the critical 
frequency f^ be defined as 



the frequency ratio u be defined as 


CO _ 2nf _ 

1^0 ~ ^ 0 ° fo 


(D5) 


(D6) 


and the frequency response excitation number as 
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Then equation (D2) can be written as 


ot. 




P( 0 )) = 


(“o + 


2r Qf^ 

* 0 

77^ + (U - 1)^ 


L ’I +1 


7)^ + (U + 1)^ ’ 

7)^ + + 1 J 


This function is comrnonly plotted logarithmically. Let 

PSDrgai(f) = 10 lo^io “o> *^0 ’ 

where 


(D7) 


(D 8 ) 


(D9) 


Cq ^lOlogjo 


(DlO) 


E(u, 7 ),c) = 10 logjQ [E^(u, 7 j,c)Eg(u, 77 ,c)] 


(Dll) 


E^^UjTJ.c) 


7)^ + (u - 1)^' 
. 77^+1 . 


(D12) 


^g(^7 C) 


T}^ + {\X+ 1)^ 
+ 1 _ 


(D13) 


The first two terms of equation (D9) determine the level of the PSD function at zero fre- 
quency ratio and are independent of frequency. The last term will be referred to as the 
basic frequency response function. The behavior of E(u, ? 7 ^c) is shown as a function of 
frequency ratio u for various excitation values r} in figure 8 where -E(u, 7 ?, c) is 
plotted. 

Figure 8 indicates that for n less than 1. 0 the value of t] determines the sharp- 
ness of the peak of -E(u, t), c) with smaller values of r] producing higher and more nar- 
row peaked curves. As the value of r) increases to one the peaks become broad and 
lower. At high frequencies curves with j} less than 1 merge to form a single curve. 
Curves with 17 greater than 1 have parallel slopes at high frequencies. 
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APPENDIX E 


BASIC SPECTRAL DENSITY FUNCTION 


The particular realization of the spectral density function defined in appendix C will 
be formulated from the reciprocal of the Fourier transform of a "damped cosine oscil- 
lation" derived in appendix D. The resulting formula is 




^’real^St) = 


1 = 


a^ cr(f2,a.) 


N. 


a(n) 

b(J2)| 


12 


\=1 


(El) 


where < N^, the complex numbers a^ and bj^ have negative real values and posi- 
tive or zero imaginary values, and a^ is a normalization constant. 

From equations (D5) to (D13) the PSD function can be written in terms of nondimen- 
sional parameters u and t] as 


PSD 


real 


(St) = Cq + 


M N. 

^ E(u,7),aj) - y E(u,77,b^) 
i=l k=l 


where now the critical Strouhal number is defined as 


(E2) 


St ^ 

The quantity u(c) is the Strouhal number ratio defined as 


(E3) 


u(c) = — ?L_ 

Sto(c) 

The quantity rj is the Strouhal response excitation constant defined as 


(E4) 


r}{c) = 


Re{c) 

^m{c) 


(E5) 
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and 


Co = 10 logjo 



(E6) 


Thus equation (El) can be interpreted in terms of the sum of a level deter minir^ 
constant and various curves such as are shown in figure 9. 
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APPENDIX F 


SIMPLIFIED SPECTRAL DENSITY FUNCTION 

Equation (El) is the power spectral density of a real random process. A modifica- 
tion of equation (6), which yields a simpler expression that is equally suitable for curve 
fitting PSD'*’ data will be discussed next. (However, although the new equation can be 
used to satisfactorily fit a spectral density, it does not represent a power spectral den- 
sity of a real random process since it is not an even function of frequency. ) 

The modification consists of dropping the term (Q - c* )/((? + ycc* ) from an equa- 
tion similar to equation (D3) of appendix D in which jto is replaced by Q, (as done in 
appendix C). The equation for the calculated lossless spectral density then 

can be written as 




I 

^0] 

i 

a 

p (f2 - a^^) 
=1 

k= 

' (n - )\) 


a.{a) 


b(f2)| 


(FI) 


where all the symbols are the same as in equation (El). 

Equation (FI) also can be written in terms of nondimensional parameters as the sum 
of terms. The equation for the calculated lossless PSD function is then 

a b 

^^^cal^^*^ " ^0 2 F(u,77,a.) - ^ F(u,? 7 ,bj^) (F2) 

i=l k=l 


where u, 77, and are defined by equations (E 3 ) to (E6) and 

F(u,r/, c) = 10 log^Q E^(u,t), c) (F3) 

where E^(u,7j,c) is defined by equation (D12). 

In figure 9 the basic and simplified Strouhal response functions are compared for 
Strouhal response excitation 7 ? values of 0. 05 and 5. 0. The figure indicates that the 
two functions have very similar shapes. 

Equation (El) can be made into an even function of frequency by using the absolute 
value of instead of Thus 
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M, 


12 


S| I ^1^1 ’ V 

i=l 


Nx 


TTd^^l- V 

k-:. 


(F4) 


Equation (E4) represents a physically realizable (real random) process. However, 
since no use of the function at negative frequencies will be made herein, this notation 
where lfl( replaces O will be dropped. 
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APPENDIX G 


TRENDS OF RATIONAL FUNCTION AND CURVE FIT PARAMETERS 


The following will be discussed: the effect of velocity on PSD^^^; the variation of 
PSD^g^j with microphone angle; the effect of the effective exhaust velocity on OASPL^; 
the variation of OASPLp with microphone angle. 

Effect of velocity on PSD^^^. - A series of PSD^^^j^ curves that pertain to meas- 
urements made at 2. 74 and 5. 48 meters (9 and 18 ft) at the 100° angle over an effective 
velocity range of 162. 8 to 242. 3 meters per second (534 to 795 ft/sec) are shown in fig- 
ure 10. These PSD^^I^ curves are based on data at both microphone heights. The data 
at the microphone height of 1.37 meters (4. 5 ft) were not used in these calculations 
since they contain irregular variations not shown in the data at 2. 74 or 5. 48 meters 
(9 or 18 ft; see fig. 5). Also shown in figure 10 are the corresponding PSD'*' points. 

The functional form used for the PSD^^j^ curves is 


PSD+^l . 10 logiQ. 


{Cl - b^)(n - b2)(J2 - bg) 


(Gl) 


Figure 10 shows that equation (Gl) fits the data well. The deviation of the data from 
the empirical curve is less than 1 decibel at most frequencies. The maximum of each 
curve occurs between Strouhal numbers of 0. 25 and 0. 35. A problem was encountered 
in determining the maximum, since the initial Strouhal number for the lower velocity 
cases is very close to that at which the maximum for the PSD^^^ curve occurs. In 
figure 11 the curves shown in figure 10 are plotted together and fall within a narrow 
band. 

The variation of PSD^^^ with angle. - A series of curves which pertain to 

measurements made at the maximum effective velocity (V^ = 242. 3 m/sec (795 ft/sec)) 
over a range of angles from 20° to 120° is shown in figure 12. 

Equation (Gl) was used to curve fit the data. Since only data from the microphones 
at a height of 2.74 meters (9 ft) was available at most angles, the curves shown in fig- 
ure 12 are based on less data than those in figure 10. 

Most of the PSD^^^ curves shown in figure 12 are similar. For instance, they all 
tend to peak at a Strouhal number between 0. 25 and 0. 35. 

Two importance characteristics of the PSD^^^ curve are the Strouhal number at 
which it peaks and the peak value. The variation of these two characteristics with 
microphone angle is shown in figure 13. Here, in order to indicate that the trends 
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applied to more than one particular case, results based on data from the two highest 
velocity tests are shown. Note that the peak of the Strouhal curve is calculated from 
one -third -octave SPL data; thus it can represent only the approximate maximum. Also 
note that the peak is not as well defined for the 224 . 6-meter-per-second (737 -ft/ sec) 
data as for the 242. 3 -meter -per -second (795-ft/sec) data since the number of points 
with Strouhal number smaller than that of the peak is less. 

In figure 14 the calculated curves for angles from 40° to 120° of figure 12 

are plotted together. This figure shows that all the curves fall within a narrow band. 

Effect of jet velocity on OASPL^. - Figure 15 shows the variation of OASPLp 

with effective velocity at each angle. The same trends observed in reference 18 are 
present. The OASPL"^ from the microphones between 20° and 100° can be fitted by 
lines having a sixth power slope. The OASPLp from the microphones at 120° show a 
slope closer to eighth power. This is especially apparent at the highest velocities. 

These results suggest that the spectrum from the microphones between 20 and 
100° is dominated by blown- flap noise, since flap noise is thought to be proportional to 
the sixth power of the effective velocity and that the microphones at 120° are affected by 
jet noise, since jet noise is proportional to the effective velocity to the eighth power. 

Tabulated results. - To enable the reader to use the actual empirical curves, the 
coefficients found are listed in table H. Also in the table are values of the cost function 
Crj,, the OASPL"^ value, and the reflected- ray transfer function parameters Q and 5. 
The ranges of each of the parameters seem limited, but they were not analyzed for 
trends since most of the data used to obtain them contained little information on the low 
frequency roll off. 
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TABLE I. - EXPERIMENTAL CONDITIONS FOR TEST OF WING 
FLAP SYSTEM BLOWN BY TURBOFAN ENGINE 
WITH COANNULAR NOZZLE AND TAKEOFF 
FLAP SETTINGS OF 0^, 20°, 40° 

[Horizontal position of wing, 30. 48 cm (12 in. ); vertical position of 
wing, 104. 14 cm (41 in.); fan area, 5096,8 cm^ (790 in.^); core 
area, 1812,9 cm^ (281 diameter, 93.8 cm (3.077 ft); dis- 

tance from source to microphone, 30. 48 m (100 ft); microphone 
height,^ 2.74 m (9 ft).] 


Core 

Fan 

Effective 

Ambient 

Humidity, 

velocity, 

velocity, 

velocity, 

temperature 

% 

V 

core 

Vfan 

V 





m/sec 

ft/sec 

m/sec 

ft /sec 

m/sec 

ft/sec 

°C 

°F ' 


263.7 

865 

231.6 

760 

242.3 

795 

11. 67 

53 

54.9 

236.2 

775 

219.5 

720 

224.6 

737 

12.22 

54 

53.2 

216.4 

710 

207.3 

680 

210.0 

689 

12.22 

54 

53.2 

193.5 

635 

192.0 

630 

192.3 

631 

11.67 

53 

57.6 

157.0 

515 

164.6 

540 

162.8 

534 

11.67 

53 

54.9 


^At 100° microphones were also at heights 1/2 and 2 
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TABLE n. - PARAMETERS CALCULATED FOR EMPIRICAL EQUATION {Gl} 


Effective 

Total 

a 


^1 

b 

1 





Overall 

Magnitude of 

6 

velocity, 

cost 








sound 

reflected-ray 
















function. 


Om 


u 



Qe 

9m 

ke 

9m 

pressure 

transfer 



m/eec 

^tot 











level, 

function, 














OASPL 

Q 



Microphone angle, 0^ 



242.3 

0. 159 

419 

375 

-1.37 

1. 06 

-252 

191 

-2.13 

1.29 

-0.69 

1.59 

108.2 

0.206 

8.4 

224.6 

.179 

465 

375 

-1.20 

.98 

-248 

195 

-2-14 

1.47 

-.712 

1.5 

106.2 

.051 

7. 

3 

210 

.22 

563 

356 

-1.09 

.975 

227 

259 

-2.0 

1-65 

-.693 

1.48 

104.8 

.002 

5. 

98 

192.3 

.31 

480 

567 

-.431 

.944 

-198 

285 

-2.04 

1,3 

-.592 

1.5 

102.8 

.13 

6. 

28 

162.8 

.36 

644 

538 

-.0815 

.96 

-199 

265 

-2.28 

1.14 

-.0836 

1.14 

99,63 

.117 

6.28 







Microphone angle, 20 

o 






242.3 

0. 0215 

469 

— 

371 

-1.18 

1.14 

-254 

895 

-2.37 

1.33 

'0.748 

1.5 

111.7 

0.518 

0. 

574 

224.6 

.018 

456 

375 

-1.28 

1.14 

-255 

113 

-2.26 

1.39 

.78 

1.51 

109. 7 

,483 

5. 

72 

210 

.036 

465 

375 

-1.24 

1, 07 

-256 

109 

2.3 

1.38 

-.718 

1.39 

108. 1 

,425 

5.65 

192.3 

.368 

548 

403 

-1.097 

1.30 

-186.4 

266.3 

-1.97 

.832 

.531 

1.55 

106.82 

.118 

6. 

28 

162.8 

.144 

457 

517 

-.7 

.91 

-238 

154 

-2.35 

1.63 

-.578 

1.33 

102. 1 

.378 

5. 
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Microphone angle, 40 

o 







242.3 

0.374 

418 

378 

-1.35 

1.17 

-262 

-629 

-2.34 

1-59 

-0.77 

1.51 

113.7 

0.489 

5. 

75 

224.6 

.0388 

411 

379 

-1.46 

1.33 

-256 

-72. 8 

-2.32 

1.53 

-.746 

1.47 

111.6 

.405 

5. 

71 

210 

.043 

449 

375 

-1.36 

1. 13 

-253 

-95 

-2.37 

1.63 

-.726 

1.39 

109.4 

.398 

5. 

69 

192.3 

.047 

421 

378 

-1. 46 

1.033 

-263.6 

-79.3 

-2.38 

1.8 

-.752 

1.43 

107.6 

.401 

5, 

75 

162.8 

,057 

428 

372 

-1.53 

.890 

-265 

-92.9 

-2.39 

2.3 

-.78 

1.42 

103.7 

.417 

5. 
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Microphone angle, 60^ 


! 

242.3 

0. 019 

448 

375 

-1.24 

1.041 

-261.2 

-29, 86 

-2.39 

1.33 

-0. 758 

1.50 

114.7 

0. 483 

5, 

66 

224.6 

,0249 

436 

389 

-1.26 

1.2 

-250 

-55 

-2,39 

1.03 

-.72 

1.47 

112.3 

.352 

5. 

82 

210 

.048 

459 

375 

-1.34 

1. 19 

-254 

-62 

-2,36 

1.28 

-.698 

1.39 

110.44 

.372 

5. 

73 

192.3 

,0536 

499 

375 

-1.13 

1.33 

-225 

-116 

-2.63 

.874 

-.512 

1.46 

108.1 

.315 

6. 

28 

162.8 

.07 

512 

345 

-1.42 

1,35 

-233 

-111.3 

-2,49 

.988 

-.213 

1.27 

104.1 

.205 

6, 

28 



Microphone angle, 70^ , 

242.3 

0. 034 

598 

365 

-0. 81 

1 

-248 

-328 

-2.95 

0.57 

-0.525 

1.54 

113, 8 

0.239 

6. 

28 

224.6 

.0365 

562 

331 

-.96 

1.15 

-234 

-96.8 

-2.79 

.646 

-.552 

1.61 

111.8 

.254 



210 

.0477 

520 

345 

-1.05 

1.09 

-241 

-91.1 

-2.55 

.843 

-.58 

1.42 

110. 46 

.198 



192.3 

, 054 

503.7 

350,3 

-1.23 

1..07 

-245 

-100 

-2.4 

1.126 

-.606 

1.33 

108. 18 

.148 



162.8 

.072 

500 

345 

-1.4 

1.12 

-244.4 

- 109. 1 

-2.42 

1.45 

-.314 

1.22 

104. 48 

. 06138 

_ 






TABLE II. - Concluded. PARAMETERS CALCULATED FOR EMPIRICAL EQUATION (Gl) 



Effective 

velocity, 

m/sec 

Total 

cost 

function, 

^tot 





b 

1 

b 

2 

*^3 

Overall 

Magnitude of 
reflected-ray 
transfer 
function, 

Q 


6 

He. 


He 

Urn 

He 

9m 

He 

9m 

He 

9m 

sound 

pressure 

level, 

OASPL 

Microphone angle, 80® 



242.3 

0. 0356 

586 

559 

-0.538 

0. 896 

-241 

-95.1 

-3.35 

0.597 

-0.512 

1.45 

114.25 

0.259 

6. 

28 

224.6 

.0411 

620 

339 

-.691 

.955 

-231.2 

-113 

-2.96 

.558 

-.483 

1.51 

112.36 

.224 



210 

.055 

596 

344 

-.714 

1.05 

-226 

-110 

-2. 88 

. 491 

-.475 

1.52 

110.4 

.244 



192.3 

.073 

548 

354 

-.978 

.961 

-244 

-95.9 

-2.81 

1. 06 

-.523 

1.31 

108.5 

.143 



162.8 

.082 

558 

352 

-1.21 

1.02 

-246.2 

-84.8 

-2.52 

1.26 

-.188 

1.14 

104.7 

,123 






Microphone angle, 90® 



242.3 

0.008 

499.7 

374 

-0.71 

0.674 

-258 

672 

-2.5 

2.36 

0.83 

1.28 

113.8 

0.302 

5. 

73 

224.6 

.0123 

476 

370 

-.827 

.68 

-264 

9.8 

-2.5 

2.26 

-.828 

1.3 

111.57 

.287 

5. 

71 

210 

.0145 

499.6 

374 

-.94 

.78 

-258 

-33.6 

-2.57 

2.17 

-.764 

1.28 

109.4 

.269 

5. 

47 

192.3 

.0472 

516 

370 

-.92 

.84 

-253 

-48.68 

-2.73 

1.95 

-.613 

1.26 

107 

.0988 

6.28 

162.8 

.056 

501 

368 

-1.09 

.814 

-254. 5 

-68.9 

-2.53 

2.09 

-.535 

1.22 

102. 86 

.0964 

6,28 



Microphone angle, 100® 


242.3 

0. 00872 

472 

380 

-0. 768 

0.694 

-262 

106 

-2.31 

2.36 

-0.809 

1.34 

114.2 

0. 0997 

5. 

94 

224.6 

.0123 

457 

375 

-.732 

.674 

-259 

59 

-2.39 

2.51 

-.832 

.1.29 

111.2 

,139 

6,14 

210 

.0138 

491 

374 

-.909 

.782 

-253 

99.6 

-2.47 

2.62 

-.703 

1. 18 

108,9 , 

.0731 

6, 

21 

192.3 

.01 

503 

376 

-.938. 

.747 

259 

51.86 

-2.56 

2.39 

-.736 

1,21 

106.1 

,102 

5. 

63 

162.8 

.017 

509.3 

373.2 

-.99 

.782 

-257 

-11.2 

-2.64 

2.706 

-.534 

1.15 

102. 05 

. 108 

5. 

67 


Microphone angle, 120® 


242.3 

0.0158 

368 

375 

-1.42 

0.748 

-288 

113 

-2.17 

— 

2.26 

-0.839 

1.58 

— 

115.6 

— 
0. 07 

7. 

66 

224.6 

.013 

396 



-1.197 

.647 

-279 

111. 2 

-2.2 

2.46 . 

-.839 

1.48 

112.3 

.058 

7.28 

210 

.014 

443 



-1.04 

.685 

-269 

114.2 

-2. 19 

2.29 

.827 

1.39 

109.6 

,0399 

7. 

93 

192.3 

.0127 

437 



-1.19 

.743 

-270 

66.6 

-2.28 

2.47 

-.784 

1.23 

106.4 

.0439 

7. 

77 

162.8 

.0346 

460.2 

374.5 

-1.15 

.735 

-265 

33.6 

-2.48 

2.76 

-.709 

1.21 

101.42 

.0614 

8. 57 

Microphone angle, 140® 

— 

242.3 

0.0274 

382 

375 

-1.88 

1.21 

-288 

-79.9 

-2.25 

2.66 

-0.727 

1.24 

114.7 

0. 0102 

7. 

2 

224, 6 

.034 

432 

375 

-1.56 

1.028 

-279 

-18.6 

-2.35 

2.2 

-.76 

1. 23 

m 

.0276 

5. 

18 

210 

.039 

464 

328 

-1.53 

.896 

-277 

106 

-2.22 

2.25 

-.714 

1,22 

110 

.029 

9, 

3 

192.3 

.043 

383 

375 

-1.91 

.791 

-287,5 

-111.2 

-2.38 

3.84 

-.59 

1. 35 

108. 98 

,130 

5. 

12 

162,8 

.047 

417 

375 

-1.76 

.711 

-281 

-89.96 

-2.45 

4. 08 

-.529 

1.31 

104. 47 

.056 

4. 628 


cn 

CO 






ib) Turfaofan engine under*the-wing in externally blown -flap test configuration at takeoff flap setting of 0®-20’^-40^, 



<c) Microphone array. At 100^ measurements also taken at 1/2 and 2I1 q. 
Figure 1. -Test facility geometry (ref. 21. 



Sound pressure level, dB 



Frequency, Hz 


(c I Microphone angle, 100°; no apparent cancellations or (d) Microphone angle, 120°; no apparent cancellations or 
reinforcements. reintorcements. 

Figure 2 . - Typical one -third-octave sound power level spectra from under-th e-wing, turbofan engine, externally blown-flap 
test configuration (ref. 1). Effective exhaust velocity, 242,3 meters per second (795 ft/ sec); distance from source, 

30.48 meters (100 ft); microphone height, 2.7 meters (9 ft). 
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Strouhal normalized pressure spectral density, PSDjg|. dB Calculated unaffected SPL^gi, dB 



(a) Comparison of calculated and measured 
sound pressure level. 



id) Calculated Strouhal normalized pressure spectral (el Comparison of PSD^ based on measured data corrected 

with the hypothetical reflection effect (PSOj.) and PSD'^ 
as calculated from curve fit (PSD^g|). 


Figured. - Initial and final iteration values of SPL, SPL^ |T(f )|2 and PSDcai. Microphone 6(]Ofrom inlet. Effective 
exhaust velocity, 242. 3 meters per second (795 ft/ sec); distance, 30. 4S meters (100 ft); microphone height 
2.74 meters 19 ft); engine centerline height, 2.74 meters (9 ft). 
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Sound pressure level, d8 



50 100 200 500 1000 2000 5000 10 000 20 000 


(a) Microphone height, 2.74 meters <9 ftl, 



Frequency, Hz 

(b) Microphone height, \.V meters (4.5m. (ci Microphone height, 5.49 meters 118 ftl. 

Figure 5. - Comparison of SPLj-^i, SPL^p and SPL^g| one-third-octave spectra at three microphone heights. Microphone 
angie, 100°, Effective exhaust velocity, 242.3 meters per second (795 ft/ sec); distance, 30.48 meters (100 ft). 



Frequency, Hz 

Figure 6. ' Comparison of SPL(,g|, SPL^f,. and 
SPL^gl for iargest magnitude of reflected-ray 
impedance ( Q =0.51 evaluated by computer 
procedure. Microphone angle, 20°; effective 
exhaust velocity, 242. 3 meters per second 
(795 ft/ sec); microphone distance, 30.48 me- 
ters UOOft); microphone height, 2,74 meters 
(9 ft). 
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Magnitude of reflected-ray impedance. Q, dimensionless 








Figure 9. - Comparison of basic and simplified frequency (Strouhal) response function at frequency (Strouhal) 
response excitation (t?) value of 0.05 and 5.0. 
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.2 .4 .6 .8 1 2 4 6 8 .2 .4 .6 .8 1 2 4 

Stroulial niimber, St, fOg/Vg 

tb) Effective exhaust velocity, 224. 6 meters per second (737 ft/ seel. (c) Effective exhaust velocity, 210, 0 meters per second 

(689ftfsec). 

Figure lO. - Comparison of caiculated with measured PSD"^ for different effective exhaust velocities. Microphone angle, lOQP. 









Strouhal number, St, fOg/V^ ' 

*"‘9ure a - Effect of velocity on calculated PSO Jg, for microphone 













{al Variation of peak Strouhal tb) Variation of peak value, 

number. 


Figure 13. - Variation of peak Strouhal number and peak value Of calculated 
PSD'*’ with microphone angle. 


cH 



Figure 14. - Variation of calculated lossless 
spectral density. Effective exhaust velocity, 
2A2 . 3 meters per second (795 ft/ seclj 
turbofan engine externally blown-flap 
configuration; flap setting, 0°-20®-40®. 
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Effective exhaust velocity, Vg, m/sec 

(g) Micrc^hone <h) Microphone it) Microphone 
angle, 100^. angle, 120®. angle, 140®. 


Figure 15. - Variation of lossless overall sound pressure level 
with effective velocity. Turbofan engine externally blown- 
flap configuration; flap setting, 0®-20®-40®; microphone 
radius, 30.48 meters (100 ft). 
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